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Abstract 

Retaining only the 'timelike' component A of the vector potential a 
skelet model with explicit global center symmetry is constructed for SU(2) 
Yang-Mills theory. It is shown that the A gluon vacuum is equivalent with 
the 4-dimensional Coulomb-gas. In 1-loop approximation, the effective 
theory of the skelet model exhibits non-local self-interaction. A coupling 
constant A is found that governs the loop expansion. For A < 1 the ef- 
fective theory does not confine (the string tension vanishes), whereas for 
A = 1 confinement (with non-vanishing string tension) takes place. In the 
deconfined phase the A Q gluons exhibit non-zero rest mass and global center 
symmetry is broken by the vacuum state. Confinement sets on for A = 1 
when the rest mass of A gluons vanishes and the global center symmetry 
of the vacuum is restored. 
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I. INTRODUCTION 



Lattice QCD results indicate that confinement and global center symmetry are strongly 
correlated |], @, |3|, ||J]. Simulating pure Yang-Mills systems on the lattice it was found: 
(i) there exists a critical temperature T c for quarkless QCD at which phase transition 
takes place; (ii) the free energy of a static qq pair rises linearly with increasing separa- 
tion distance and the string tension obtained from the Wilson loop and that from the 
correlation of Polyakov lines agree well for T < 0.6T C ; (in) the confining phase main- 
tains global center symmetry, but it is broken in the deconfined phase. It is the most 
important message of the lattice results for us that the transition from the confined 
to the deconfined phase seems to be accompanied by breaking of the global center 
symmetry. In the present article we investigate the continuum SU(2) gauge theory 
with global center symmetry. We show that the confined phase (with non-vanishing 
string tension) may exist with global center symmetry, whereas in the deconfined phase 
(with vanishing string tension) this symmetry is broken. In connection with the phase 
structure of the theory we also show that only a few of the infinitely many parameters 
are of significance. 

Our treatment can be outlined as follows. We start with the generating functional 
of the Green's functions of the lattice regularized SU (2) Yang-Mills theory. In order to 
ensure global center symmetry the ST/ (2) invariant Haar measure is used in the path 
integral [Q, and the gauge fixing is performed without violating global center symmetry. 
This leads to a cut-off dependent periodic effective potential. For the sake of simplicity 
we neglect the spatial components of the vector potential remaining invariant under 
global center transformations. The 'skelet' SU(2) Yang-Mills theory obtained in this 
way contains only the timelike Aq components of the vector potential being sensitive 
to the global center transformations. 

We show that the skelet model is equivalent with the macrocanonical ensemble 
for a, D = 4 gas consisting of infinitely heavy point charges interacting via Coulomb 
interaction. This structure of the vacuum is completely induced by global center sym- 
metry. In addition to the strong coupling constant the model contains infinitely many 
other parameters: the fugacities of the charges ±1, ±2, etc. 

We expand the effective potential at its minima and derive the effective action 
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in 1-loop approximation ||. The 1-loop order effective action contains a non-local 
self-interaction. A parameter A is found governing the loop expansion. 

In the framework of the effective theory we derive an expression for the string 
tension. The string tension is determined through the propagator of the gauge field 
in the effective theory. According to the value of the loop expansion parameter A the 
propagator of the effective theory has qualitatively different behaviour. For A < 1 it 
exhibits simple poles whereas for A = 1 a double pole occurs. In the first case the 
string tension vanishes, while it has a non-vanishing value in the second case. 



II. SKELET MODEL 



The generating functional Z[j] of the SU(2) Yang-Mills theory is by definition the 
vacuum to vacuum transition amplitude in the presence of the external source j": 

Z\3\ ~ jvA-J VAq exp j-ij / d'xF; u F; v -S gf - I d*xj«A; - S Bh } (2.1) 

with the gauge fixing term S g f and the ghost term S g h- We modify the generating 
functional of the usual perturbative treatment in order to ensure its invariance under 
global center transformations following @, f|, Hf, but we take global center transfor- 
mations into account in a rather explicit way and use gauge conditions which do not 
violate global center symmetry. 



1. Under an arbitrary finite SU(2) gauge transformation h(x) the spatial and the 
timelike components of the vector potential transform as follows^, [|: 

(A a y) = h(x)(A^h-\x) + h(x)d i h- 1 (x), 
(expjia^y})" = (expjia^y})/^). (2.2) 

The center of a group consists of the elements commuting with any element of 

the group. The center of the group SU(2) consists of the elements Zq = 1 and 

Z\ = —1. Global center transformation means the same center transformation at 

every point of 3-space at a given time slice. This transformation takes advantage 
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of the fact that we can perform finite gauge transformations at any time slice 
separately. In the continuum limit this leads to field configurations discontinuous 
in time. The global center transformation Z\ = — 1 = exp j^vr 1 ^} is not trivial. 
Parametrizing the general group element h(an) = exp{ia;n|} by the colour 
vector a = an (a G [0, 2ir], and n arbitrary unit vector in the solid angle An), we 
get from Eq. ( |2.2j ) for the effect of the center transformation: 

h[an)z x = h (-{2tt - a)n) . (2.3) 

The colour vector an flips to — (2tt — a)n due to the non-trivial center transfor- 
mation Z\. In order to treat the field configurations obtained by non-trivial global 
center transformation explicitly, we allow for global center transformations Zi at 
each time slice t (k t = for identity transformation and kf — 1 for z\ = —1). 
Taking the parameter of the local gauge transformation from the sphere of 
radius 2tt, agt £ SW, and taking the sum over k t = 0, 1 we cover the SU (2) group 
twice. The multiple counting of the gauge equivalent intermediate vacuum states 
shall be removed by dividing the path integral by the appropriate multiple of the 
volume of the gauge group. 

The physical vacuum to vacuum transition amplitude Z[0] is constructed in 
Hamiltonian formalism starting with the gauge Aq = (a = 1,2,3) 0, || (see 
Appendix A). The whole interval T of time evolution is divided into small inter- 
vals of length a. At each time slice the projection to colour singlet vacuum states 
is performed by the operator: 



iz kt ' 



J2 I V H (hz k <) (hz k >) = £ fv H h{hz 

!=0,1 fct=0,l 

Ell/ d H^t h ((agt + 27rk t (-l) ht )n gt ) , (2.4) 



fe t =0,l x 



d H a = d 3 a— ^w. (2.5) 



with the gauge invariant Haar measure: 

1 sin 2 (±a) 
4ir 2 a 2 

(d 3 a the flat measure). Finally the path integral is divided by the multiple of the 
volume of the gauge group: 

Ell/ d n3m = [ V u h. (2.6) 

{fc i= 0,l} xt JS ^ Jg 



Later we make use of the integral 




with arbitrary functional $[aA 0a ]. We can reintroduce the zeroth component of 
the vector potential M via 

A« it = ga% t /T. (2.8) 

2. We introduce the gauge fixing conditions F a [AQ a ,yl^ a ] in a gauge invariant way 
by inserting 
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1 = A[A] / Vh]\5[F a [A h ]] (2.9) 

■* 5 a=l 

in the path integral. The functional A[A] is invariant under gauge transforma- 
tions including global center transformations as well. Making use of the gauge 
(and center) invariance of the integration measure, and the A functional, we can 
factorize out the multiple of the group volume jg Dh and can divide with it. Thus 
we find for the physical vacuum to vacuum transition amplitude Eq. (|A.7| ): 

Z[0] ~ / V(aA a ) f VA1 AL4]| Fs0 (f[ 5[F a [A]]) JW*™*^ (2.10) 

* ' \a=l ) 

with the action Q]: 

S[A h ( A °^,Ai) = -—J^Jd'xF^F^, (2.11) 
and with the generalized field strength tensor: 

Ht a = -\ (4Sa-4?)> (2.i2) 

For infinitesimal gauge transformations the generalized field strength tensor takes 
the usual form. 

3. We complete our definition of the vacuum to vacuum transition amplitude by a 
particular center invariant choice of the gauge conditions. According to (|2.3|) the 
colour vector Aq does not change its orientation in colour space due to center 
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transformation. Further on colour vectors aAnn and (aA — 2n)n are center 
equivalent. Therefore the gauge conditions 

F'lA] = aAl m = 0, F 2 [A] = aA 2 0St = 0, 

F 3 [A] = F{aAl)$ t+a - F(aJ§) B t = (2.13) 

do not violate global center symmetry for J 7 (a) periodic, J 7 (a) = JF(a + 27r). The 
gauge condition F 3 = is the discretized version of the condition ad J 7 (aA ) = 0, 
and can be considered as the center invariant generalization of the temporal 
gauge OqAq = 0. This is a generalization allowing for finite jumps of the zeroth 
component of the vector potential at any of the time slices due to global center 
transformation. As for small fields we would like to recover the temporal gauge, 
we require J 7 (a) ~ a for a — > 0. The gauge condition F 3 = exhibits as a 
rule many solutions, many Gribov copies exist. To avoid this complication we 
require that the Gribov copies are those field configurations obtained by global 
center transformation. It is achieved if the equation J 7 ((3) — J 7 (a) = has only 
the single solution (3 = a for a, (3 6 [0,27r]. For example the function J 7 (a) = 
tan (hex) satisfies all these requirements. In the derivation of the vacuum to 
vacuum transition amplitude, however, we have to make use of the properties of 
the function J 7 (a), but we do not have to specify it. The sums over the global 
center equivalent configurations (over kts) are just those over the Gribov copies 
in the path integral. 
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The A[A] functional defined above takes the following explicit form (see Appendix 



A[A]\ F ^ = exp|TV r ln(47T 2 )-^ilnsin 2 QaAo^+X:in|^ / Mo 



xt/ 



xt 



(aA, 



xt 



(2n-aA z t ) 2 



(2.14) 



da 



J 7 {a). Inserting this into the path integral 2[0], we carry out the integration 



over Al and A^ and also the sums over fc t 's (see Appendix C). Making use of the global 
center invariance of the function J 7 (a) and that of the action, we can replace S(F 3 ) by 
(see Eq. O ): 

2 sin 2 [\aA 3 : 



^ (Gxt+a - Qxt) 
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(2.15) 



with G{a) = | (a — sin a). This allows us to recover the Haar measure and eliminate 
the dependence on the function jF(a). Shifting the argument of the Dirac delta by 
a field independent value cgt and integrating over it by Gaussian weights, we obtain 
(Appendix C): 

Z[0] = fVuf VA a t exp I iS[A h{u) j , A j ] + ^ In 2 + ^ In sin 2 (~au) 

1 /l 1 \ 2 1 

where au = aA^ G [— 2n, 2n]. 

We see that retaining global center symmetry has lead to the periodic tree level 
effective potential due to the Haar measure, as expected^, |3|, ffl. 



We see from Eq. ( |2.2| ) that the spatial components of the vector potential re- 
main unaltered under global center transformations. For our goal is to investigate the 
consequences of the global center symmetry we neglect the spatial components of the 
vector potential by setting A® = in the present work. 

In this way we obtained the 'skelet' model which only contains those degrees of 
freedom being sensitive to the global center transformations. In Euclidean space the 
generating functional of the Green's functions takes the form: 

Z[q] = Z[q]/Z[% (2.17) 

with 

Z\q\ = J Vue- Ss - s *f-f dixgu , (2.18) 
and a 3 jg = q the external source coupled to u, where 

^ j4 /oi„.\2 



S s = ^Jd 4 x(d l uY + V[u}, (2.19) 

S B f = ^-Jd*x±(d g(au)) 2 , (2.20) 
and the tree level effective potential is given by: 

V = J d^xViuix)), V(u) = --j In sin 2 (-auj - In 2. (2.21) 



As the action S s + S g f exhibits u — > —if symmetry, we can rewrite the integrals: 

r2n/a _ fin/a rriT ^ _ 4i _ f2n/a 

t 

l-2n/a 



/2ir/a r2iv/a . /'27r/a 

du x ...^\{2 du x ... = e TV «-^ 2 U du x ... (2.22) 
-2w/a ~ JO ~ JO 



The additional exponential factor just cancels the zero frequency term of the Fourier 
expansion of the tree level effective potential. 

As the tree level effective potential has minima at au = ±7r, the field configura- 
tions close to these contribute most significantly to the path integral. Therefore we use 
the expansion 



g(a) = c(±)r) + * 



{a T tt) + 0((a =F vr) 3 ) 



Q = ±7T 

1 



= aT^ + 0((aT7r) 3 ), (2.23) 



^ : 2^ 



which choosing £ = 1 leads to: 

^ | rf 4 x(9 oM ) 2 . (2.24) 
Then we obtain for the generating functional: 

= J Vue- s[u] 'I d4xqu , (2.25) 

with the action 

S[u] = ^Jd 4 x [{d l uf + (d u) 2 ] + VoM (2.26) 

and with the tree level effective potential 

V = V-TVa- A \n2 = J d 4 xV (u). (2.27) 
The Fourier expansion of the periodic effective potential is given by: 

V Q (u) = Vv CT cos (craw), (2.28) 
a ^o 

with v a = V- c = | cr | 1 . The tree level effective potential is periodic in the interval 
— 2n < au < 2-k with minima at au = ±7r, and it tends to +oo for au = 0, ±2n. 
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III. ANALOGY WITH THE COULOMB GAS 



Having performed the Fourier expansion of the effective potential and the Taylor ex- 
pansion of exp{— Vo}, we carry out the Gaussian path integral over the field u and 
obtain[||: 

Zc[q] = E { —T- II E E Mi) zc[q], (3.1) 



n=0 \j=l Xj a (j)¥=0 



with 



zc[q] = (DetD) 1/2 exp | ^ E |E if7 0')^ - q x J a 2 D xy - 

• (x>(*)<W"fc)}- (3-2) 
For vanishing external source 

Zc[0] = (Detr 1/2 Et^fflE E v aU \ex V {-W°U)Dx 3 x k <r{k) 

n=0 U - \j=l Xj a (j)^0 J \ Xj*k 

(3.3) 

is identical with the macrocanonical partition sum of the D = 4 dimensional Coulomb 
gas with infinitely heavy (static) point charges a — ±1, ±2, .... The Fourier coefficients 
— v a of the periodic tree level effective potential play the role of the fugacities of the 
charges o. The free u-propagator, 

a 2 g 2 

a 2 D(xj,x k ) — > — — for a — > (3.4) 

n 2 (Xj - x k ) z 

represents the Coulomb interaction between the static charges. In addition to the 
strong coupling constant g the model contains infinitely many other parameters, the 
fugacities v a . 

We conclude that the requirement of global center symmetry leads to the equiva- 
lence of the A gluon vacuum with aD = 4 dimensional Coulomb gas of static charges. 
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IV. EFFECTIVE THEORY 



For now on we use lattice regularization which is rather natural from the point of view 
of the definition of the path integral and of the dependence of the effective potential 
on the lattice size a. Although the calculations have been done in lattice regularized 
form, for the sake of simplicity we write the expressions in their continuous form if it 
is not misleading. The main steps of obtaining the effective theory are as follows: 

1. The effective potential at the tree level has minima at au = ±7r. Therefore the 
vacuum expectation value (u(x)) = u does not vanish as a rule, but corresponds 
to one of the minima of the effective potential. Let us shift the field variable 

u(x) = u + T)(x), (4.1) 
and write the generating functional of the Green's functions as: 
Z[q\ = e ^- A j^x q J Vr]exp l J_ J d 4 xr]n v 

-a' 4 J d 4 xr]q - J d 4 xV (u + rj) J . (4.2) 

2. Now we can linearize the exponential of the path integral in terms of r)(x) by 
introducing the auxiliary field h(x) by the definition: 

exp / d 4 ^n??j = (DetD) 1/2 • 

■J X>/iexp|^ J d 4 xd 4 yh(x)a 2 D(x,y)h(y) + J ( ^h(x)ar](x)^j . 

(4.3) 

The generating functional has now the form: 

Z[q] = (DetD) 1/2 e a " a_4 / d4 ^y > P/ie^^ d4 ^ 4 ^ Wa2 ' D(:E ' y)h(y) • 

• J Vr]e a ~ 4 f d4xr i( h -q)-f d 4 xV{u+r)) ^ ^ 

In order to clarify the meaning of the auxiliary field let us go back to Minkowski 

space, which roughly means the replacements d 4 x — > —id A x, r\ — > irj, q — > iq, 
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h — > i/i, and forget about lattice regularization. Then the naive integration over 
T>rj leads to the Dirac-delta functional: 
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h(x) — q(x) + ^2 cr(j)a5(x — 



x 



(4.5) 



According to that the auxiliary field h(x) can be interpreted as the effective 
charge density being the sum of the external charge density and the polarisation 
charge density. The latter arises due to the Coulomb charges in the vacuum. 



3. We expand the tree level effective potential in the Fourier series ( |2.28| ) and the 
exponential e~ v in Taylor series, and perform the path integral over rj. In order 
to integrate over all group elements in terms of the new variable r] x , we choose 
the interval of integration as ai] x £ [— 7T, tt]. Then we obtain: 

Z[q] = e aUa ~ 4 J d4x i J jyfa^ f dixd4 yH x ) a2D ( x >y) h (y)+ lnz l h -q} (4 

where 

eMMd = (DetD )i/* £ t}L ( / ^2 J2 vJ\lS*[h-q], 

n=0 ,L - \j=l J a a(j)^0 ) 

(4.7) 

with 

/(»)[£] = g'EM^Jjr^e^^^ 



■ y m c • -(./)• ( 8) 

and = h x — q x . Here we used the correct discretized expressions of the lattice 
regularized theory. We can write: 

lnz[£] = ln/o[£]-m2o[£], (4.9) 

with 

in/oK] = E lnDet (° 2£, ) + E[ msiiih K-)- m ^]- ( 4 - 10 ) 
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The functional z [£] has the form of the macrocanonical partition sum of an ideal 

gas: 



n - 

oo / i \n 



-in = ln E— r EE^ 



-e 



kr(au+7r) 



n=0 n! 



& + ia 



= -ee^f4— ek7(as+ ^ ( 4 - n ) 

According to that we can interpret the Coulomb gas of charges a as an ideal gas 
of quasiparticles in the external field h x . 

4. To the next we perform the path integral over h(x) by the saddle point method 
taking into account the tree level and the 1-loop contributions. The solution of 
the saddle point equation 



_5_ fl 

Sh. 



\^Y.h x a 2 D xy h y + \nz[h-q]^ =0 (4.12) 



is a local functional h 0x [q}. Performing the saddle point integration over the 
auxiliary field h x , the generating functional can be written as: 

Z\q\ = e- M X.«- Zo [q]Z 1 [q], (4.13) 

where Z Q [q] and Zi[q] are the tree level and 1-loop contributions, respectively: 

Zo[q] = exp h 0x [q]a 2 D xy h y[q]+ In z[h [q]-q]^, (4.14) 



Z x \q\ = [Det(a 2 Q)} V \ (4.15) 
with 



a Q xy = a D xy + — . (4.16) 

OllQxOllOy 

5. The generating functional W[q] = In Z[q] of the connected Green's functions 
takes the following form: 

W[q\ = const. +W [q} + W 1 [q], (4.17) 
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where Wo [9], and W\[q\ are the tree level and 1-loop contributions, respectively, 
for the Coulomb gas, 

W [q] = -auY,Qx + ^^2h 0x [q]a 2 D xy h y[q\+ In z[h [q}- q\, (4.18) 

x xy 

W x [q] = -^TrHa 2 Q). (4.19) 

6. Let us derive the effective action in 1-loop order defined as the Legendre trans- 
form of the generating functional W[g]: 

m = -T, a (Xx + u)q x -W[ql (4.20) 

X 

where the classical field Xx is given by 

axx + au = JJpSi (4.21) 
Sq x 

as the functional of the external source q x . 

7. In order to make our treatment self-consistent we require 

(au x ) = au, i.e. a Xx\ q=0 = 0. (4.22) 

Furthermore we require that the effective potential has minimum for the constant 
field configuration u x = u. In this way we ensure that the action is expanded at 
the classical solutions. 

The solution of the saddle point equation ( |4. 1 2| ) and the u values satisfying the 
self-consistency equation ( j4.22| ) were found (Appendix D). It was established in 1- 
loop order that the effective potential has extrema at au = 0, ±7T, ±2ir. au = ±7r 
correspond to the minima of the tree level effective potential. Therefore we 
perform the expansion of the theory around these configurations and later check 
that the 1-loop effective potential has minima for them. It is a consequence of 
global center symmetry and that of the charge conjugation symmetry u — > —u, 
that the minima of the effective potential remain at the same place if quantum 
fluctuations are included. 



The solution of the saddle point equation ( |4.12|) is given by 

1 
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hox-qx = -{Q, 1 Dq)x+ l -U(Q G \Q Q l Dqf) +0(q 5 ), (4.23) 



with 

.2 



a 2 Qo) xy = a 2 D xy + f 2 5 xy . (4.24) 
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We construct the effective action T[x] by inverting the functional dependence 
given by Eq. (|4.21|) and inserting it in the r.h.s. of Eq. (}4.20|) retaining terms up 
to the fourth order in the classical field. 

Although the skelet model contains infinitely many parameters, the fugacities v a 
and the strong coupling constant g, the effective action in 1-loop order can be 
expressed in terms of the parameters 

f 2 = 2 • 2! £ v a (a)- 2 > 0, U = -2 ■ 4! £ v a {a)- 4 < 0, (4.25) 

CT>0 cr>0 

and the new coupling constant 

A = -7^>0. (4.26) 

This is a significant reduction of the number of relevant parameters of the theory. 

Calculating the 1-loop contribution Wi[q\ to the generating functional of the 
connected Green's functions in Appendix E, we show that it is of order A. There- 
fore we investigate the skelet model for < A < 1, where the loop expansion 
seems to be justified. The 1-loop contribution turned out to be convergent only 
if s = 2f2/g 2 < 1/8, so s is a small parameter in the model. The expressions 
listed below are obtained by simple algebraic manipulations from the tree and 
1-loop contributions to the generating functional, Wo[q) and W^ifg], resp., given 
in Appendix E. 

The classical field Eq. Q4.21 ) as the functional of the external source is given by: 



X, = -/ 2 ((l + A-/ 2J D- 1 + 0(A 2 , S 2 ))g 



-\hqt + i/4/ap-V). + - A hf2ql{D~ l q) x + 0(X 2 ,q 5 ). (4.27) 
oo 4 

Here f ^D -1 is of the order s. For we did not go beyond the 1-loop approximation, 
the operator in the linear term on the r.h.s. is known only up to 0(A) terms. 
Therefore we must take its inverse with the same accuracy. Then we get for the 
inverse relation 



q.i 




4/: 



xl{D-\) x + 0{X\x 5 ), (4-28) 



which was found by iteration first solving the equation without the self-interaction 
terms. The effective action takes now the form: 

m = \ E X, (^S xy + D^j Xy - £ Xi (4.29) 

Here we included only the leading order local terms of the self-interaction, taking 
into account that ~ @( s )- 

We see that f 2 > 0, and / 4 < are the correct signs for the effective potential 
having minimum at Xx = 0. That justifies the expansion made. 

The effective action describes a non-local theory with quartic self-interaction. We 
read off the inverse propagator of the A$ gluons from the quadratic term of the 
effective action: 

V-\p) = ^ + D~\p) + 0(\ 2 ,s 2 ), (4.30) 

J2 

with the lattice regularized inverse propagator of the bare theory[Q: 

f 2 D-\p) = 4 £ (1 - cos( ai f)) . (4.31) 

A rest mass is generated due to the periodic effective potential, i.e. the global cen- 
ter symmetry of the effective action. This rest mass can be read off by comparing 
the — > limits of the inverse propagator with the general form g~ 2 (p 2 + m 2 ): 



m 2 = a - 2 -{l-\) = a- 2 ^ l —^>Q. (4.32) 
s Si »A 



In the neighbourhood each of its minima the effective potential at 1-loop order 
has the following form (neglecting the off-diagonal non-local terms): 

ViM = ^V^fe) 2 -^Ete) 4 (4-33) 

It is interesting to note that the rest mass as well as the quartic self-interaction 
vanish for A = 1 in 1-loop approximation. We show in the next chapter that 
the propagator exhibits rather different analytic properties for A ^ 1 and A = 1, 
which lead to a non-vanishing string tension in the latter case. 
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V. STRING TENSION 



Let us determine the free energy of a static quark- ant iquark pair positioned on the z 
axis (the quark at z — 0, the antiquark at z — L). In the SU(2) Yang-Mills theory the 
Wilson operator is given by: 



w 



(L;[A}) = Trfe-^^A^ e -i/>| fl A S r- 

V 2 = 



<-o=T 



z=L 



x°=0 



(5-1) 



which leads in the skelet model to the expression: 



w(L;[u\) = Trexpji^ dx ^ (u(x°, 0, 0) — u(x°, 0, L)) r a 1 , (5.2) 



i.e. 



w 



(L; [ X ]) = 2 cos 2 + X*)) = 2 cos ^ £ , (5.3) 



with the source 



Qx — S x o5 y o(S Z L — 5 z q). 



(5.4) 



The free energy of the static pair is defined by the vacuum expectation value of the 
Wilson operator: 



_ HL) JV X w(L; [ X ])e- r M 

with the approximate form of the effective action: 

1 



^2Xx{D )xyXy 



xy 



The Gaussian integrals can be performed easily: 

e -HL) = 2e ~^y QxVx 



(5.5) 



(5.6) 



(5.7) 



Then we express the propagator in coordinate space through its Fourier-transform and 
take the sums over x and y. This results in the following expression for the free energy: 



d 3 p 



sin' 



--a 3 

2 a a Jb (2n) 3 V- 1 (p ± ,p z ;p° = 0) 
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(5.8) 



where the 3-momentum integral is taken over the Brillouin zone. 



Calculating the string tension we replaced the propagator T> by T>o, i.e. neglected 
the quartic self- interaction vanishing for A = 1. 

The momentum integral can be performed noticing that the integrand exhibits 
poles on the complex p z plane for vanishing transverse momentum p± = 0. Therefore 
it is a good approximation to set p± = 0: 

HL) = jLx, (5.9) 
where the remaining one-dimensional integral, 

x = ha r" dp = h r dv pM (5 ,o) 

J-w/a 1 — A + S — S COS{ap z ) S J -TV ( — cosv 

with v = ap Z) C = 1 + (1 — A)/s = 1 + m 2 a 2 /2 > 1, has significantly different 
properties for A < 1 and A = 1. 

Expressing the denominator in terms of exponentials, the momentum integral 
takes the form T = I + 1 + + I_, where 

lo = — , X± = -— / e » . 5.11 

2sJ-tt( — cosv 4s J -n ( — COS V 

The integrand has simple poles on the complex v plane at v — iv± = i ln(£ ± y/(' 2 — 1) 
(v + > 0, f_ < 0) corresponding to the zeros of the equation ( = cosru± = coshw ± . We 
calculate these integrals by closing the integration path through the lines C\. (7r + i0) — > 
(tt + iR), C 2 : (tt + iR) -> (-7T + iR), C 3 : (-tt + iR) ->• (-tt + iO) for J and J+, and 
through C[: (tt — iO) -> (vr-LR), C 2 : (tt-LR) -> (— tt — ii?), C£: (— tt — -> (-vr-iO) 
for T_, where i? — > oo. We show that the contributions of these lines to the integral 
vanish for L/a =integer, i.e. in the lattice regularized theory: 

k( f + / ) dv = k ( r _ r\ dv = o (5 12) 

2s \Jd Jc\J ( — cosv 2s \Jo Jo / ( — cosv 
and using cos(=F7r ± ia) — — cosh a, 

_£(/ + t\*L^ = Jir**-^to(h,\ = n ( 5.i3) 

4s \Jci JCi) ( — cosv 4s Jo C + cos ha \a J 
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for L/a =integer. Similarly the contributions from C[ and C 3 cancel for X_, and 

— / 7 ~ e -> 0, 

2s JC* 2 C ~ COSf 



f 2 f dv iLv L R _ R 

— / — e a ~ e a — ► U, 

4s Jc 2 C — cosf 

h 1 dV e -i*«~ e -^*-0 (5.14) 



4s Jc^ C ~ cos v 

for _R — > oo. Then we can deform the integration path into the circles C+, and C_ of 
radius p — > centered at the pole iv + for X and X+, and at ii>_ for X_, resp.: 

M — 



2s Jc + C ~ cosy' 
h I dv 



J + = -±±<k —1 e & 

4s Jc+ C — cost" 

X_ = -A/ ^ e -i#«. (5.15) 
4s ic_ C ~~ cosf 

We write w = \v± + pe 11 ^, and ( — cosv = isinhf ± pe lip + \Cp 2 e 2lLp + C(p 3 )- 
Now we have to consider the cases A < 1 and A = 1 independently. 



For A < 1 we have C, > 1 and v± ^ 0, and the nominator is of the order p. Retain- 
ing in the denominator the terms of the order p, i.e. writing exp {±i^pe llp j ~ 1, 
we get: 



Jo i si 



2 s Jo isinhf + pe lLp s sinht> + ' 

f 2 [2* pidcp l h * l 

_L + = — — / —e a + = e <• , 

4s jo i sinh v + pe lip 2s sinh v + 

f 2 f 2n pidtp l_ f 2 7T l_ 

X_ = — / T-r-. — e« = — — : — e« . (5.16) 

4s Jo i sinh i>_ pe 1 ^ 2s sinh t>_ 

The free energy of the static quark- ant iquark pair is then given by: 

HL) = T^— f-^-^T— + — i— e^V (5.17) 
47ra s \smhi> + 2smhf + 2smht>_ J 

There is an attraction between the quark and the antiquark, which tends expo- 
nentially to zero for increasing separation distance L. A kind of Debye screening 
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takes place governed by the rest mass m of the Aq gluons through v± and £. 
Consequently the string tension vanishes: 

k = lim -^ = 0. (5.18) 

For A = 1 we have ( = 1 and v± — 0, and the nominator of the integrand on 
the r.h.s. of the integrals Eq. ( [5 .15 ) becomes of the order p 2 . The simple poles 



at v — iv± merged into a single double pole at v — for A — > 1. Then we 
have to retain the terms up to the same order p 2 in the denominator, writing 
exp jii-^pe 1 ^} 1 ± i^pe lif . The terms of the order 1/p vanish, because they 
contain the integral of the function e~ 1Lp for a period. Therefore Tq = and we 
get for the free energy: 



T TT f 

F{L) = — (X + +J_) = — ^. (5.19) 
Ana 4a 2 s 



We see that the free energy exhibits the 'area law'. Now we obtain a non- vanishing 
string tension: 

-2/2 ~29 2 -2 ~fi /r on\ 

k = a — = a — = a — — , (5.20) 



4s 8 8/ : 



2 



where we used A = 1, i.e. — / 4 = f 2 g 2 in the last equation. The simple propor- 
tionality of the string tension k to the strong coupling g 2 will be modified by 
taking the quartic self-interaction as perturbation into account. 



It is interesting to note that the 1-loop effective potential ( |4.33| ) becomes constant 
(the coefficient of the quartic self-interaction term vanishes) for the limit m — > which 
means that the vacuum state does not break the global center symmetry. Oppositely, in 
the case A < 1 the effective potential has two minima and the vacuum expectation value 
(u) must be one of them, so that the vacuum state violates global center symmetry. The 
result obtained is an indication rather than a proof for the existence of the confining 
phase in the model, because the 1-loop approximation can be problematic for A ~ 1. 
Anyhow, the 'area law' of the Wilson loop is a consequence of the treatment which 
does not violate center symmetry. Therefore we see that global center symmetry is of 
extraordinary importance for the existence of a phase with non- vanishing string tension 
even in the continuum theory. 
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VI. CONCLUSIONS 



Summarizing we constructed the skelet model for SU(2) Yang-Mills theory with global 
center symmetry retaining only the component of the vector potential. Due to the 
treatment not violating global center symmetry a cut-off dependent periodic effective 
potential occurs in the skelet model. The minima of this potential are fixed by charge 
conjugation and global center symmetry. It has been shown that the Aq gluon vacuum 
is equivalent with the 4-dimensional Coulomb gas. This structure of the vacuum is 
entirely induced by global center symmetry. 

The effective action for the skelet model has been derived in 1-loop approxima- 
tion. In this approximation the effective potential contains only a finite number of 
parameters which are particular combinations of the strong coupling constant and the 
fugacities of the static charges. The parameter A is found which governs the loop ex- 
pansion. Depending on its value, the vacuum can be in either deconfined or confined 
phase. 

For A < 1 the string tension vanishes, and the renormalized rest mass of the Aq 
gluons given by Eq. ( |4 . 3 2| ) is positive. Consequently a test SU{2) charge is screened 
at a distance 1/m due to Debye-screening. The effective potential is periodic and its 
curvature at the minima is defined by the renormalized mass m. The ground state is 
characterized by the non- vanishing expectation value of the gluon field (either au = it or 
— 7r) and the global center symmetry is broken. Thus the theory decribes the deconfined 
phase for A < 1. 

It is a hint on the possibility of the existence of a phase with unbroken global 
center symmetry that the renormalized rest mass vanishes for A = 1 in 1-loop approxi- 
mation. Then the simple poles of the effective propagator of A gluons on the complex 
momentum plane merge into a single double pole and this leads to a non-vanishing 
string tension. The test charges are not screened now and the effective potential be- 
comes a constant. Therefore the global center symmetry broken for A < 1 is restored 
for A = 1 in the ground state. These are just the properties of the confined phase. 

Finally some remarks on our procedure and the results obtained above: 

• Our procedure of deriving the effective theory contains a non-perturbative resum- 
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mation of infinitely many vertices generated by the periodic effective potential, 
i.e. by global center symmetry. This results in a vacuum structure being equiva- 
lent with a 4-dimensional Coulomb gas. Due to reserving global center symmetry 
during the whole treatment we expect our approach to reveal the IR (long dis- 
tance) features of SU(2) Yang-Mills theory correctly. On the other hand, it is 
necessary to incorporate the spatial components of the gauge field for the exten- 
sion of the model to the UV regime. 

• Our gauge fixing does not violate global center symmetry and is done rather 
consistently than in our earlier treatment |7j . We have shown that the periodic 
effective potential is completely determined by the Haar measure. It is indepen- 
dent of the explicit form of the periodic gauge condition F 3 , as one expects^. ||. 

• Part of the self-interaction (higher order in the field) and the terms higher order 
in the loop expansion were neglected. The qualitative features of the effective 
potential determined by global center symmetry are not sensitive to this approx- 
imation. The critical value of the parameter A for which confinement sets on can 
be modified by the self-interaction and higher order terms in the loop expansion. 
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A Vacuum to vacuum transition amplitude 

Let us denote by H the Hamiltonian of SU(2) Yang- Mills theory in the gauge Aq = 
(a = 1,2,3). Dividing the interval T of the whole time evolution into iV intervals of 
length a and inserting the projector operator Pg of the colour neutral physical states 
at each time slice, we write the physical vacuum to vacuum transition amplitude as 
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the product of infinitesimal physical transition amplitudes: 
(0|e- idT |0) phys = 



'N-l r . / nnn, a \ \ 

n En/</ s *- ^ n(4£r 

j=l /c 4 =0,l x ^ V X J I k=0 



-iaH 



O, (A.l) 



where h(ag ti ) = ho(oigti)zi% h G SU(2). Using the explicit form of the infinitesimal 
transition amplitude ||, 



A h(a) ja 
xt k+1 



-iaH 



A ja 

^xt k 



exp j-ia— / d 3 yF^ a 



(A.2) 



t—t k 



and reintroducing the zeroth component of the vector potential by Eq. ( |2.8| ), we obtain: 
(0|e^ T |0) phys = Jv H (aA 0a )JvA^ AHAo)3 ' AJ \ (A.3) 



Now we introduce the gauge conditions F a by inserting 1 in the form of (2.S) 
Dividing by the multiple of the group volume, we get: 

S V H (aA°") J VA™A[A] J g Vh' (flti S[F a [A h '}}) e^ AHAo) W 



<0|e-^|0) phys 



SgVh' 



(A.4) 



Making use of the fact that the integration measure, the functional AL4], and the 
action are invariant under any gauge transformation h', we can write the denominator 
in the following form: 

jvh'jVn (aA h ' 0a ) J VA h '™A{A h '] (f[ <5[F a [^']]) (A.5) 

Further on we make use of the transformation law of the parameters of the gauge 
transformations (see Eq. 



h(A )ti = h(A%). (A.6) 
Then we can remove h! from the integration variables A h °, and A h> J : 

(0|e^ T |0) phys = Jv H (aA 0a ) J VA ia A[A] (f[ 5[F a [A}}) e^ AHAo) M'l. 

* \a=l ) 

(A.7) 
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B A[A] functional 



Let us write the general gauge transformation hg t e SU(2) in the form: 



(B.l) 



where af3g t G [0, 2tt], and n t = 0, 1 integer at the time slice t. (For the sake of simplicity 
we omit the lower indeces (x, t) of the variables if it is not confusing. ) 



According to Eq. (|2.2| ), the following general expressions are obtained: 



( A o°f 



iS in(ia^° 

w 

in [\aA h ° 



A h ° 



.sin(|a^° 



A h ° 



pi fl ^0 . n2 

p cos — - — hp sm ■ 



aA \ sin (§< 







R 2 aA . . aA ^ sin (± 
p cos — p sm — - 1 



2 

aA 



P 3 







aj3 Aa aA a(3 

sm h — — sm cos — 

2 A 2 2 



^, (B.2) 



and 



h\l 



«) 2 



■l) nt + 2%n t /a 
sin (ha r 



P 1 cos ( -a nt ) + p 2 sin ( -a 



A%°(-l) nt + 27rn 4 /c 



sin (|a„ t ) 



P 2 cos ( -a nt ) — p 1 sin ( -a 



(A 



/i\3 



sm 



(l« nt ) 







COS I -Q!nt I Sm 



4 3 

+(-ir7 



sm | -a nt | cos 



2 

a/3 



a/3 



sin ( 



mi) ( \a 





sin (§a/?) 







(B-3) 



with a nt = aA (—l) nt , where we made use of the equations: 



{aAt 

sin (^A/ 

/l /' 
cos 1-clAq 1 



n A a 

aA a + 2nn t (-ir—±.. 



{-l) nt sin ( -a nt , 



{-l) nt cos f -a nt 
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(B.4) 



The expressions for sin (clAq /2^J, i.e. for Aq" in terms of /5 a, s we find by taking the 
sum of the squares of the expressions given by (|B.2|) : 



sin 



2 a -^0° 



l_(l + cos 2 0) sin 2 ^ 



sin 2 Qa/?) 



2 Uj4 -o 



, - cos^-r^- sin(a/?) sin(a/lo) 

2 2 Aq 



+ sin 

= B(a(3,0), 

where the spherical coordinates 

(3 1 = (3 sin 9 cos ip, (3 2 = (3 sin 9 sin </?, /3 3 = /3 cos 6 1 
have been introduced. Then we get: 

clAq 10 = arccos[l-2S(a/3,0)]. 



(B.5) 



(B.6) 



(B.7) 



The inverse of the functional A [A] is given by: 
A-'IA] = Yl Jv H {aF) f[ 8[F a [A h °* nt \\ 



(rl-K e-K r2n \ Sin 

TT / d{a(3 3t ) (a(3g t ) 2 / sm9g t d9 St / dipg t —r — - , 

%Xh Jo Jo 4tt 2 (a/^) 2 

3 

■Y[5[F a [A hoznt }}. 



(B.8) 



a=l 

Here we can write: 



SIF^A* *?]] = J] 



xl 



V(-l)"«+2^/a sjll sin ^ sin A flA ( _ 1)nt 

iaA n h 0(-l)"0 ^ 2 ' ^ 2 



= n 



sin(±aA£°(-l)"t) 

8{(p-(p v ) 



.it 



AQ°(-l)™t+2nn t /a . (\ 



where the equation 



cos QaA (-i; 



sin 



sin# 



(B.9) 



(B.10) 
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has 1 root for ip u G [0, 2n] and n t = of the form ip u = [aA (— l) nt — (2u + l)vr]/2 with 
v = for 7r/2 < aA /2 < n, and v = — 1 for < aA /2 < n/2. For n 4 = 1 we can 
change the integration variable from (p to —ip and then we get the same root as for 
n t = (notice that | sin(aA /2 — ip u )\ = 1); 



S[F 2 [A hoz 



nE 

xt /u=l,2 



= nE 

M=l,2 



= nE 

ju=l,2 



5(0 - 0,) 



5(0 - 0,) 



A n °(-ir t+27r n t /a sin /i N og 
sm| iaA$°(-l)»«J V2 7 

5(0 - M ) 



sm 



(B.11) 



sin 

where sin^ = has 2 roots, 0^ = 0, 7r. Now we perform the integrals over the angles: 
A" 1 [A] = 



( 



= E 

{nt} 



,2- 11 sin 2 UaA$>(-l) nt ) 



i[F 3 [(/ 02 "^]] 



(B.12) 



For M = 0, vr and A\ = A 2 = 0: 

B(a[3,6^ = 0,7r) = i (1 - cosa M ) = sin 2 ^ 
with a M = a(A|] + /9cos0^), and then we obtain: 



a, 



for |a M | G [0,27r] 



and 



a^° 3 = 



sm | 



(B.13) 



sm | 



27r — for G [27r, Air] ' 



a M for G [0, 2ti] 

a M — 27rsgna M for G [27r, 47r] 



(B.14) 



(B.15) 



24 



As the function J 7 (a) is invariant under global center transformations, the sums over 
rifS can be carried out, and the inverse of the A functional takes the following form: 



rifS can 

/ 

1 1 Al=A 2 =0 Z - 



{«*} 



i-2-k 



r** 1 1 sin 2 [f (^ + /3cos^)' 

st Jo m r _i ), 



aA* 3 ) 



- F(aA%° 3 " 



^0=^0=0 




1 



) 2 + (2vr - aA h ° a y 



1 



F{aA h Q ° 3 ) 



t+a 



Inserting the expression (|B.15| ) we find: 

A-MA] , , = TT d(a(3 St ) -^-Y sin 2 (-(Al + (3 cos 

1 1 

a(Ag + /3cos0 M )| 2 [2tt- \a(A 3 + (3 cos 



(B.16) 



0, 



A' 



•5! 



a(A 3 + (3 cos # M )| 2 [2tt- |a(A3 + /5cos^)|] 2 . 
J^(a(A 3 + /3 cos ^)) | t+a - H<Al + P c os 0J>) Q . 

(B.17) 



We introduce the new integration variable: £ = aft cos 8^. Then the integrand 
becomes independent of /i and the sum over /i can be carried out: 



E 



n cos Bp ju 



2n cos t 



2tt 



-2tt 



(B.18) 



Furthermore we have to use the field configurations satisfying ^{aA^)^— F^aAfyt = 0, 
i.e. ay4|| t = aAo| t+a , and av4Q| t+a — 27rsgn y4Q| t+a . As JF is periodic, the Dirac delta 
can be rewritten in both cases in the same way: 



5(F(aA 3 + Ot +a -HaAl + a 



F(aA 3 t + £ t 



(B.19) 
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According to the requirements satisfied by J- {a), the argument of the Dirac delta has 
2 zeros in the interval £ [ _ 27r, 2ti\. ^ kt = £ t+a , £ t+a — 27rsgn£ t+a . For both zeros hold 
the equations 

1 



sin 2 Q( a ^ t + ^ t )) 
F\aAl t + tikt) 



sin 



3 

t+a 



F\aAl t+a + t t+a ) 



(B.20) 



In the case of the second zero we can make the integral transformation from £ t+a to 
£ t+a — 27rsgn£ t+a , under which the rest of the integrand remains unaltered. Then we 
can take for both zeros: 

1 1 



Mo t+a + Sfct 



+ 



(2ir-\aA s ot+a + i 



kt\ 



+ 



(B.21) 



\aAl t+a + i t+a \ 2 (2tt- \aAl t+a + i t+a \f 

Finally both zeros give the same contribution to the integral cancelling the factor \ 
before the sum over /i. We integrate over the successively according to their time 
argument t simply writing £ kt = £ t+a : 



A" 1 [A] 



F=0 



= n 



xt 



1 sin 2 ^{aA* gt + & T/2 
4^ \r(aA 0£t + ^ T/2 )\ 

1 

+ 



aA l st + 1 (2vr - \aAl St + ^ /2 |) J _ 
Choosing £jt/2 - * for the pure gauge final state, we obtain: 

1 sin 2 (|a^ t ) / i i 



(B.22) 



A- 1 [A] 



F=0 



n 

xt 



+ 



An 2 \r(aA 0St )\ {(aA^,) 2 (27r-|a^ a |) 



(B.23) 



C Explicit form of the vacuum to vacuum transi- 
tion amplitude 

We have made the jumps of the vector potential A explicit in the path integral. This 

resulted in the sum over center equivalent copies of any field configuration. As the Haar 

26 



measure, the A functional, the action, and the gauge condition F 3 are invariant under 
global center transformations, only the arguments of the Dirac deltas with the gauge 
conditions F 1 and F 2 transform in a non-trivial way under global center transformation. 
Taking this into account, we obtain: 



Z[0] 



= E / v H{aA a Q ) [ VA m A[A] F = (f[5 F a 

{n t } JS ^ J \o=l L 



. AS[A h o( A o) i 



F=0 



= e y ^tM) y ^ 

{n t } 



27m f 



^xpjE.lnsin 2 QaAo) -E.ln(47T 2 )} 
rLMo + 27m t (-l)™0 2 



A[A] F= „ 



>3A3] 



aAn 



.6(F(aAl) t+a -F(aAl) t )]e iS l Aho{AoU 

v /" t) t„A*\ r v v exp ln sin2 fe^) ~ Ex ln(47r2) 1 

• exp {E ln(4vr 2 ) - £ ln sin 2 QaA ) + ]T ln |^(aAg) 1 1 • 



+ 



_(aA ) 2 (27r-aA)) 2 
■ (jl8(aAl) 5 (aAl) 5 (^(a^) t+a - ^4*) t )) ^^'M'] 

We perform now the sums over n t 's: 

*[°] = /% a t^y^(n 5 (^Mokt + a-^Mok) 

. e iS[A h o^o) 3 , A j] + J2 x ln |^'(aAg)| | 



F=0 



(C.l) 



IF=0 



(C.2) 



In order to recover the effective potential induced by the Haar measure, we rewrite 
the Dirac delta on the r.h.s. and show that the physical vacuum to vacuum amplitude 
is independent of the choice of the gauge function J 7 (a). The argument of the Dirac 
delta exhibits 2 zeros: aA^ t = aAl t+a , and aAl t+a — 2-K^^a.A\ t . Consequently, we can 
write: 



5(F(aAl) t+a -F(aA 3 ) 
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= \r(aA 3 0t )\ \- 6 1 ~ aA ° t+a ) + 6 ' ~ ^ t+a + 2nsgnA ° t+a)] ■ ( C - 3 ) 

We make the integral transformation t+a — > Aq t+a — 27rsgnAQ t+a in the second term. 
Under this transformation the rest of the integrand remains unaltered as J 7 , J 7 ', and 
the action S remain invariant. Then both terms give the same contribution to the path 
integral and we can replace the Dirac delta as follows: 

5 (HaA 3 ) t+a - H«Al) t ) -> * 6(aA* t -aA* t+a ). (C.4) 

Let us take now the function Q{a) = \{ol — sin a), being monotonically increasing as its 
first derivative G'{ot) = sin 2 f > 0. Then the equation Q(f3) = G(a) has only a single 
solution p = a. Making use of this property, we can write: 

2 sin 2 (ha A 



6 (F(aAl) t+a - F(aAl) t ) -> j^ajjjp * (^Mot) ~ GMtJ) ■ (C5) 

Inserting this into the physical vacuum to vacuum transition amplitude, the terms 
Z) x ln|.F / | cancel in the exponent of the integrand, and we recover the Haar measure 
induced tree level effective potential: 

*[°] = / ^flat^o / ^ (U S ^( aA lht + a - «) a )j • 



•e 



Shifting the argument of the Dirac delta by a (vector potential independent) 
constant c x and integrate over c x by Gaussian weights, we obtain: 

/ Vc exp j-^ ^ c 2 1 . . . (j{ 6 (Q{aAl t+a ) - Q{aA\ = 

= exp|-^i:^(^Mo t+ J-^ t )) 2 }--- (C7) 
Then the physical vacuum to vacuum transition amplitude takes the form: 
Z[Q] = Jv m J%JvA$^™W + ^^.*>^(W). 

•^P {-^ £ ^ (S(a^ t+ J - S(a^)) 2 } . (C.8) 

The additional interaction term expressed through the function Q{a) is also defined 
completely by the Haar measure. 
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D Solution of the self-consistency and the saddle 
point equations 



We find the solution iteratively in two steps. At first we solve equations ( 4.12|) and 
( 4.22|) in the tree approximation. Then we show that the functional ho x arid the values 
u are not modified by the 1-loop corrections. 



1. In tree approximation the self-consistency equation (|4.22| ) 

6\nz[g\ 



(D.l) 



£x=hox[q=0] 



is consistent with the saddle point equation ( |4.12| ), 



0. 



(D.2) 



£x=ha x [q]— q x 



only if the functional /ioxlV] vanishes for vanishing external source, i.e. h 0x [q 
0] = 0. In light of that we expand ln,z[£] in Taylor series at £ x = 0: 



OO 1 

lnz[£] = MOl+EErjM*, 

x k=l K - 



(D.3) 



where fk is the k-th functional derivative of In z with respect to £ x at £ x = 0. We 
take into account the terms up to the order 0(C, 5 ). It is easy to show that the 
odd derivatives of ln/ vanish at £ x = 0, 5 2k+1 (lnl )/5^ x k+1 \^= = 0. From Eq. 
( P . 1| ) we get fi = 0, which then means 



olnz [£] 



cr^O 



1(7 



ia(au+n) 



(6, + Kr)" 
cos(awcr) + isin(attcr) 



(D.4) 



<T^0 



Since f CT = t>_ CT , the real part of the sum is automatically zero. The imaginary 
part, 



(-1) CT 
v a sin (aw) = 
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(D.5) 



vanishes for u satisfying sin(awcr) = for any integer a, i.e. for au = vtt, 
[y = integer). Then we get e mua = (— iy a , and 

fk = (-l) fe+1 E^(-l)^ +1) 7T^I' ( D - 6 ) 



cr^O 



i.e. all the odd derivatives of In z vanish at £ = 0, f2k+i — 0, and the even 
derivatives are given by: 

f 2k = (-1) W 2E^(-1)^ +1) ^. (D.7) 

<7>0 17 

Note that v = ±1 correspond to the minima of the effective potential in the tree 
approximation. For v = ±1 holds: 

f 2k = (-1) W -2E^, (D.8) 

cr>0 ° 

i.e. f 2 > 0, / 4 < 0, etc. We shall see that these are the correct signs ensuring 
that the effective potential at 1-loop order has minima at au = ±7r. 

We rewrite the saddle point equation ( j4.12|) : 



YF D *v + MX = -Ya 2 D xyqy - l -Ml (D.9) 

y v 

including terms up to the order (9(£ 4 ). Introducing the matrix 

(a 2 Q ) xy = a 2 D xy + f 2 S xy , (D.10) 

we can write: 

& = -Y,l( a2 Q°yV D y^-luY,l( a2 Qor 1 U& (d.ii) 

yz ^ z 

We solve this equation by iteration: 

£x = -Y^( a2 Qo)~%y a2D yzQz 

yz 

3/ V / 

+0(g 5 ). (D.12) 
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(Note that the matrices D xy and (Qo) X y are symmetric.) For the sake of simplicity 
we introduce the notations (Aq) x = J2 y A xy q y , (qAq) = J2xyQxA xy q y and do 
not write out the powers of a (they can be introduced at any step to obtain 
dimensionless expressions), then: 



-{Q^Dq) x + i/ 4 (Qo \Qo l Dqf) x . 



(D.13) 



Including terms up to the order 0(q 4 ), we obtain: 



£ = (Qo Dq) x - ^Ja(Qo Dq) x (Q (Q Dqy) 



S.t' 



C = {Q^Dq)i 



(D.14) 



for later use. 



2. Now we show that the functional ho x [q] is not modified by 1-loop corrections. 
Including the 1-loop term the self-consistency equation ( |4.22| ) is given by: 



MM + Ite'M + ^Trln(a 2 Q) 
6 2dq x 



0. 



(D.15) 



9=0 



Inserting the explicit form of the 1-loop contribution of the generating functional 
(Appendix E), 



Wife] = -\^Ha 2 Q) = ^y^Qo l Dq)l 

2 2/ 2 y 

and writing £ x [q] = ^°'[q] + h£ ^[q], we obtain for the term of order h: 



(D.16) 



q=0 



J__5_ 



Trln(a 2 Q) 



<?=o 



]__5_ 



W 1 [q] 



g=0 



k 

ft 



^ 2 ((g - 1 J D) 2 g) ; 



q=0 



0. 



(D.17) 



This means that the initial value ho x [q = 0] = is not modified due to 1-loop 
corrections. Then the solution of the saddle point equation remains unaltered 
and also the value of au is not modified. 
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E 1— loop contribution 



Now we calculate the 1-loop contribution W\[q\ to the generating functional of the 
skelet model. 



The 1-loop part of the generating functional of the connected Green's functions 
can be rewritten: 

Wx[q\ = const. - ^Trln(a 2 Q) = const. - ^Trln [{a 2 D) (l + f 2 {0 D ~ l 



const. -^Tr In (l + / 2 (0^ _1 



2 

1 oo 



const. H — 



Tr 



(E.l) 



where /2(C) is defined by ^ ^ In- 2 ^] = 8xyf2{£x)- Making use of the Taylor expansion 



(|D.3|) with vanishing odd derivatives, we get for the trace 
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+^(n - i)fr 2 (D- i r^f!(eD-^ + o(e: 



(E.2) 



and for IVifg]: 
Wife] = 



1 f 00 

const. + I fE(- 1 ) n / 2 n E^[P" 1 ) n 

4 /2 n= i x 



00 



+-^E(-i) B #E£[( £rl ) B 



+T7/ 4 2 E(-i)"t™-i)r 2 E^(^ 1 U 2 p- 1 ), p- 1 

n=l 



-l\n-2 



(E.3) 



Let us make use of the fact that the diagonal matrix elements of the powers of the 
inverse propagator do not depend on the Euclidean coordinates due to translation 
invariance (in the infinite volume limit) and reexpress the sums over n in terms of \x 2 - 



1 00 

& = 5 E(-imMi 

Z 71 = 1 



(E.4) 
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The last term on the r.h.s. of W\ we rewrite as follows: 

1 CO 

i-ji {e(D-re) E(-i)> - m- 2 ^- 1 

10 n=l 

i oo 

= rr-Jl (t\D^?e) ( L)£< T™f? ' in 

iD m=0 

i oo 

= YQ fi {e(D-re) (-i) e (-irm/rMp- 1 ) 



16 
1 

16 
1 

16 
1 

16 
1 

16 



\n-2l 



l\m— ll 



m=l 



4/ 4 2 (^(W)(-l)4 



d_ 

2 



-tEHr'r'pT']. 

m=l 

/ 2 E(- 1 ) n / 2 "[p- 1 )1 uu 

n=0 

oo 

/ 2 E(- l ) ri /2 n [p- 1 ru + /2 



(E.5) 



Here we can write: 



5 //|\ 2/ 2 



df 2 W 



(E.6) 



Finally the 1-loop contribution to the generating functional of the connected Green's 
functions takes the form: 



WiK] = const. + \f 2 ^e x + Y A f 2 ^ 
+ ^4 2 Q + 2/i 2 )E^[p- 1 ) 2 ]^ 2 . 



(E.7) 



Let us express this functional through the external source q x : 
W 1 [q] 



2/: 

+^ 2 E(^ 1 ^ 

^ 4 J 2 x 

+l9lUQo 1 Dq)l(D-%(Q^Dq)l, 
where g\ = /| (| + 2/x 2 ). 
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(E.8) 



Now we see that the 1-loop contribution is determined by /j,2" 

d^p 



/'2 



1 oo 

^E(-!) n / 2 v 

Z n=l 



b (2tt) z 



d 4 p 



a 2 D~\p) 
d 4 p 



b (2tt 4 ) 
1 



£ (1 - cos(a^)) 

At=0 



(E.9) 



Ib (2tt) 4 1 + 2f 2 g-i El =0 (1 - cos(apM)) J " 

The geometric series converges for 8s = 8 ■ 2f 2 g~ 2 < 1, i.e. s < |. We write the 
nominator of the integrand as: 4s (l — [3 X^=o cos(ap M )) with = s(l + s)" 1 . It is well 
known || that the momentum integral converges for < /3 < \ which is satisfied for 



s < 



8 ' 



Let us perform the integral over p° at first. The integrand develops simple poles 
on the complex p° plane if p —>■ |J. We estimate the momentum integral by the pole 
contribution, and set p = 0. Then the 3- momentum integration can be carried out 
trivially, and the 1-dimensional integral 



I-J2 



1 

2 



;1+4s)/3 2tt 



7r I a | 

dp - -. — -r 

-7r/o C — cos(ap u . 



(E.10) 



is left over, £ = (1 — 3/?)/3 _1 = + 1. It has simple poles at ap° = iv± = 
i In ± -\/C 2 — l) • We can take t>± ~ ±ln(2£) considering £ ^> 1 for s < ~. The 
1-dimensional integral is of the same type as the integral Xq in Part V. With a similar 
treatment, we obtain: 



/'2 



1 



1 



2 \s sinh v + 

The parameter s <C 1 is a small one in the theory. 



X 



-s. 



(E.ll) 



The generating functional W[q] of the connected Green's functions is the sum of 
the tree contribution W [q] =\nZ + const, and that of the 1-loop contribution Wi[g]: 

W[q] = -auY^qx + ^iqQiQtDq) 

x A 

(iQo'Dq) 4 ) + (?(1 - hQtQ*){QtDq)' 

1 



+ lg/4 EiQo'mKD-'tyiQo'Dq) 
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(E.12) 



with 



Q = D + f 2 = D(l + f 2 D- r ), 
Q^l + XQ^D, Q 2 = l + 2\Qv 1 D, (E.13) 



and the parameter 



\ f 4 f 4 f A /"C 1 A\ 



The parameter A governs the loop expansion of the generating functional W[q]. The 
terms proportional to A in Qi and Q 2 represent the 1-loop contributions. The last term 
on the r.h.s. of Eq. ( E. 12|) is of the order f\f 2 2 s 2 ~ f 2 X 2 and was neglected. (We 
neglected also the term proportional to /g.) 
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